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PREFACE 


Mathematical  studies  of  biological  processes  have  many  nonlinear 
functional  equations .  This  study  shows  that  the  method  of  successive 
approximations,  coupled  with  the  numerical  inversion  of  Laplace 
transforms,  provides  an  effective  computational  approach, 
y/  For  illustrative  purposes  a  nonlinear  differential  equation,  a 

nonlinear  differential-difference  equation  and  a  nonlinear  diffusion 
equation  are  considered,  a 
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It  is  shown  that  Laplace  transform  can  he  effective  in  the 


numerical  solution  of  nonlinear  functional 


t 
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I.  INTRODUCTION 

(1-3) 

In  the  previous  Memoranda  in  this  series,'  '  ve  have  demon¬ 
strated  the  effectiveness  of  the  laple-ce  transform  in  obtaining  the 
numerical  solution  of  linear  functional  equations.  In  I,  ve  con¬ 
sidered  the  renewal  equation,  in  II,  differential-difference- equa¬ 
tions,  and  in  m,  the  heat  equation,. 

Here  ve  wish  to  consider  the  possibility  of  using  numerical 
inversion  of  the  Laplace  transform  plus  successive  approximations 
to  treat  nonlinear  functional  equations.  To  Illustrate  our  tech¬ 
niques,  ve  consider  the  nonlinear  differential  equation 

u7  ~  -  u  -  u2,  u(0)  *  c,  (l) 

the  nonlinear  differential-difference  equation 

u7(t)  =  -  u(t-l)  -  u2(t),  u(t)  s  c,  •  1  s  t  s  0,  (2) 


and  the  nonlinear  heat  equation 

k(x)ut  =  \g.  +  *u2 

u(x,0)  -  c  sin  ax,  u(0,t)  =  u(l,t)  =  0  _ 


(3) 


Numerical  results  are  given,  as  veil  as  times  of  execution 
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U.  DIFHSHENl'iAL  EQUATIONS 


IPo  illustrate  the  general  approach  in  its  simplest  form,  let 
us  take  the  nonlinear  differential  Eq.  (l).  Taking  Laplace  trans¬ 
forms,  we  have 


Un')  =  -  L(u)  -  L{u  ), 


(*> 


(s+l)  L(u)  =  c  -  L(u), 
w..%  c  -  L(u2) 

*J(U)  “  ^i+1)  * 


(5) 


We  solve  this  equation  numerically  "by  mannn  of  successive 
approximations  plus  numerical  quadratures  in  the  following  fashion. 
The  function  uQ(t)  is  taken  as  the  solution  of  the  equation 

L(u)  =rfr  •  <6> 


Ignoring  the  fact  that  we  can  find  u(t)  explicitly,  we  use  the 
quadrature  techniques  presented  in  I,  H,  and  HI  to  obtain  the 
values  {uQ(ri)}. 

Next  we  write 

c  -  L(u2) 
s  +  1  * 


(7) 


To  evaluate  L(u~),  we  employ  the  same  quadrature  formula 
o 


N 

-I  vi  uo(ri)ri'1  •  (®) 

We  now  determine  using  quadrature  techniques  and  the  explicit 


inversion  formula. 
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Continuing  in  this  fashion,  ve  generate  a  sequence  of  functions 
{u  (t)]  evaluated  at  tbs  points  =  -  log  r^  i=l,2, ...,N,  vith 


Naturally,  the  procedure  will  converge  only  if  jcj  is  suffi¬ 
ciently  small.  Subsequently,  we  shall  discuss  what  techniques  can 
“be  used  in  the  general  case. 

IH.  NUMERICAL  RESULTS 

Ten  successive  approximations  for  each  value  of  c  can  he 
accomplished  in  6  sec  of  execution  time  on  the  1M  7090.  Graphs 
of  the  solutions  are  given  in  Figs.  1  and  2. 

IV.  DIFgEREmAlr-DIFgERENCi;  EQUATIONS 
Similarly,  starting  vith  Eq.  (2),  ve  have 

L(u'(t))  =  -  L(u(t-1))  -  L(u2(t)),  (10) 

whence 

1 

L(u)  (s+es)  *  c  -  J  ce"st  dt  -  L(u).  (11) 

o 


Hence,  ve  compute  uq  by  means  of 
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We  see  that  despite  the  great  increase  in  complexity  of  a 
differential- difference  equation  over  an  ordinary  differential  equa¬ 
tion,  the  computational  time  and  effort  using  this  procedure  is 
virtually  the  same  for  "both  types  of  equations.  Once  again,  ve 
have  convergence  of  this  iteration  procedure  only  if  |c|  is  suf¬ 
ficiently  small. 


V.  NUMERICAL  RESULTS 

Fifteen  seconds  of  execution  time  are  required  for  ten  successive 
approximations.  Graphs  of  the  numerical  results  are  given  in 
Fig.  3.  For  some  reasons  -which  ve  have  not  as  yet  understood,  all 
of  the  values  obtained  following  the  procedure  given  above  are  in 
excellent  agreement  with  the  actual  solution,  except  for  the  very 
last  value  of  t.  This  will  be  investigated  at  a  later  time. 

VI.  HEAT  EQUATION 

Consider  next  the  nonlinear  heat  equation 
k(x)ut  *  +  b  g(u), 

u(0,t)  =  u(l, t )  =  0,  (14) 

u(x,0)  =  h(x). 

Using  the  Laplace  transform,  ve  readily  obtain 

L(u)m  -  s  k(x)  L(u)  =  -  k(x)h(x)  -  b  L(g(u)),  (15) 

with  the  tvo -point  boundary  condition 
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L(u)  =0,  x  =  0,1*  (l6) 

Here  '  denotes  d./dx. 

As  above,  ve  employ  successive  approximations.  Let  uQ  be  the 
solution  of 

L(uo)"  -  s  k(x)L(uQ)  =  -  k(x)h(x), 

(17) 

L(uq)  =0,  x  -  0,1. 

First  we  solve  the  differential  equation  numerically  for  a  set  of 
s-values,  s  =  1,2, ...,N,  and  then  ve  use  the  Laplace  inversion 
technique.  Storing  the  values  {u  (x,^)}  for  x  =  0,A,2A,...,M&=1, 
ve  then  solve  the  equation 

1(1^)"  -  s  kfxjL^)  =  -  k(x)h(x)  -  b  L(g(uQ)). 

(18) 

b^)  =0,  X  =  0,1. 

Continuing  in  this  way,  ve  obtain  a  sequence  of  values  (un(x,r^)}. 

VH.  KU>ERICAL  RESUU5 

Five  successive  approximations  require  6  min  10  sec  of  IBM  7090 
execution  time,  and  the  numerical  agreement  was  quite  satisfactory. 
Results  of  this  method  are  shown  in  Figs.  4  and  5* 
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VHI.  C^VEtlGENCE  ASPECTS 

From  the  standpoint  of  rapidity  of  convergence,  we  would  prefer 
to  vise  quas i linearization  procedures  (cf.  Refs.  4-6).  For  example, 
in  treating  the  differential  equation  of  Section  II,  we  would 
prefer  to  use  successive  approximations  in  the  form 

L(U^)  =  -  L(un)  -  L(^.x  +  VVVi11'  (19) 


This,  however,  complicates  the  matrix  inversion  aspects.  Consequently, 
we  feel  that  it  is  better  to  use  the  simpler,  but  less  efficient 
approximation  technique. 

There  is  also  the  problem  that  the  initial  function,  the  solution 
of  the  linear  equation,  may  have  too  large  a  norm  for  convergence. 
Provided  that  the  solution  of  the  nonlinear  problem  exists,  we  can 
often  use  extrapolation  techniques  (cf.  Ref.  7). 

For  example,  we  can  replace 

^  +  u2  (20) 


by 


k<xK  =  V 


+  bu 


(21) 


■nfl  use  the  new  procedures  for  small  b.  The  solution  u(x,t)  =  u(x,t,b) 
can  then  be  determined  for  larger  values  of  b  V  using  extrapolation 


formulas 


7 


0  .4  .8  1.2  1.6  2.0  2.4  2.8  3.2  3.6 

t 


Fig,  1 


Q 


10 


4Ut*Utt+u*/IO 
u(x,0)  -  tin  T  x 
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Fifth  lttrotion 
4Ut  =  U„+u*/2 


n  n  n  o  n  n  n  o  ooo  o  o  n  n  n  on 
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*  XEQ 

*  LIST 

*  LABEL 
CNONLIN 

NCN  LINEAR  DIFFERENTIAL  EQUATION 
DIMENSION  RC15) »  W  ( 1 5 ) » A  ( 1 5  » 1 5  ) »T(15) » U ( 1 5 ) * B ( 1 5 ) 

10  FORMAT (1 12) 

11  FORMAT { 1H1 ,4X11 HD I MENS  ION  =»I3) 

13  FORMAT (//5X32HROOTS  /(5X6E20.8)) 

15  FORMAT ( //5X32HWEIGHTS  '(5X6E20.8)) 

16  FORMAT (//5X16HEXP LICIT  INVERSE//) 

17  FORMAT ( 5X6E20 * 8 ) 

101  FORMAT (//5X32HT  =  -LOG  R  /(5X6E20.8)) 

102  FORMAT ( E12 • 8 ) 

103  FORMAT (//10X4HC  =  .F5.2) 

32  FORMAT (//5X32HLAPLACE  TRANSFORM  4  /(5X6E20.8)) 

21  FORMAT (//5X32H INITIAL  APPROXIMATION  OF  U ( I )  /(5X6E20-8)) 

42  FORMAT  ( //5X23HVALUE  OF  U(I)  AT  STAGE  »  12  /(5a6E2u.8)) 

I  INPUT 

READ  10,  NA 
READ  10, N 
PRINT  11, N 
DO  12  1  =  1 »  N 
CALL  DBLRED ( X ) 

12  R(I)=X 

PRINT  13,  (R(I),I=1»N) 

DO  14  1=1 »N 
CALL  DBLRED(X) 

14  W( I )=X 

PRINT  15,  CW(I){tI  =  l,N) 

CALL  EXPINV  ( N , R  » W »  A ) 

PRINT  16 
DO  18  J.  =  1,N 

18  PRINT  17,  ( A( I , J) , J=1 ,N) 

DO  19  1=1, N 

19  T ( I ) =-LOGF { R  C I )  ) 

PRINT  101 •  (T(I),I=1,N) 

10G  READ  102  *  C 
PRINT  103,  C 

II  INITIAL  APPROXIMATION 

DO  20  1=1, N 
2u  U(  I  )  ■  C*R(  I ) 

PRINT  21,  C U C I ) ,I  =  1»N) 

III  TRANSFORM 


D 


DO  30  K=1 »NA 
DO  31  1=1 *N 
S=  I 
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3 ( I )=C/(S+1.  ) 

DO  31  J=  1 » N 

31  S ( I ) =3 ( I ) - ( W ( J ) *R ( J ) ** ( I -1 )*U(J)**2)/ (S+l  . ) 
PRINT  32 »  ( B ( I ) ,1  =  1 ,N ) 

IV  NE.V  APPROXIMATION 


DO  41  1  =  1,  N 
U( I ) =0.0 
DO  41  J=1 ,N 

41  U( I )=U( I )+A( I ,J)*B{ J) 
PRINT  42,  K> (U( I) ,1=1, NJ 
30  CONTINUE 
GO  TO  100 
END 


